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Abstract 

We extend the technique originally proposed by Honda 
et al. [1] to measure the temperature of Ytterbium and 
alkine-earth atoms confined in a Magneto-Optical Trap 
(MOT). The method is based on the analysis of excita- 
tion spectra obtained by probing the 1 So — » 3 Pi inter- 
combination line. Thanks to a careful analysis and mod- 
cling of the effects caused by the MOT light on the probe 
transition we overcome the resolution and precision lim- 
its encountered in previous works [TJ [5]. Ground state 
light shift and Rabi broadening are measured and success- 
fully compared with calculated values. This knowledge 
allows us to properly extract the Doppler contribution to 
the linewidth, thus obtaining a reliable measurement of 
the cloud temperature. We finally show how spectroscopy 
on free-falling atoms provides an alternative method to 
determine the sample temperature which resembles the 
standard time-of-flight technique. 
PACS: 37.10.De, 37.10.Gh, 37.10.Vz, 32.10.-f 

1 Introduction 

Since the first experimental realization of magneto-optical 
cooling and trapping [3], cold and ultracold atomic gases 
have become an ubiquitous tool in physics. Solid state 
physics and many-body physics [I], metrology and quan- 
tum information [5] are just a few examples of research 
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fields boosted by this new resource. Due to their favor- 
able spectroscopic properties, alkali atoms have been con- 
sidered as the most natural choice for experiments with 
ultracold vapours and they have played a prominent role 
in the spectacular advances which lead to the first observa- 
tion of Bose- Einstein condensation in atomic gases [6j [TJ [8] . 
More recently, Ytterbium and two-electrons atoms in gen- 
eral have raised interest due to their possible application 
in metrology [9] and quantum information |10) . 

Table 1: Characteristics of the 1 So — > 1 -Pi (trap) and 
1 S'o — > 3 Pi (probe) transitions. The Doppler tempera- 
ture is defined as T<j = KT /Iks- The MOT capture range 
is v c = T/k, with k = 2ir/\. According to the conven- 
tion used in the text all the quantities relative to the trap 
(probe) line are indicated by the subscript t (p). 
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Ytterbium atoms are characterized by a rich electronic 
structure. The scheme of figure [1] shows that the ground 
state 1 S'o can be optically coupled to the 1 P\ and 3 Pi 
states. Largely different linewidths are associated with 
the two transitions (see table [!} . This aspect is reflected 
in the different cooling regimes which can be achieved in 



the experiment. In the case of the 1 So — > 1 P\ transition 
a relatively large capture velocity is associated with high 
Doppler temperature. On the other hand for the 1 5'o — > 
Pi transition a much lower temperature can be achieved 
at the expense of a narrower capture range. Cooling and 
trapping of Yb atoms has been successfully demonstrated 
using both transitions Q] El El El El EES] • Combining 
these techniques with evaporative cooling in optical dipole 
traps, quantum degeneracy has been reached [TDl [TTl [T5] . 

Nevertheless it has been found that for alkaline-earth 
and Ytterbium atoms trapped by the So — > 1 Pi transition 
the measured cloud temperature is significantly higher 
than the value predicted by standard Doppler theory. Sev- 
eral explanations have been proposed to justify this dis- 
crepancy. Inelastic collisions [19] and spatial intensity fluc- 
tuations [2D] have been suggested as possible sources of 
extra heating. More recently a 3D anlysis of the MOT dy- 
namics has been carried out |21j . The predicted tempera- 
ture is in good agreement with the experimental data for 
88 Sr J22[ [23] when I/I s < 1. This approach seems to be a 
promising way to take into account the "non-Doppler" be- 
havior seen in the experiments. However it is important to 
compare the model with data obtained from experiments 
performed with other alkaline-earth atoms. 




Figure 1: (Colors online). Scheme of the lowest energy 
states of Ytterbium. The 1 S ^ 1 P 1 ( X S -> 3 Pi) transi- 
tion is indicated by a blue (green) arrow. 

In this context we extend the method originally pro- 
posed by Honda et al. Q] to measure the cloud tem- 
perature in a non-destructive way by probing the weak 
1 Sq — » 3 Pi intercombination transition. As stressed by 
Loftus et al. [2J this technique has great potential since 
it allows online monitoring of trap properties like the 



atomic cloud temperature and size, or the average mag- 
netic field at the MOT position. In this previous work the 
large linewidth of the probe laser combined with the Rabi 
broadening induced by the trapping beam limited the au- 
thors to setting an upper bound for the temperature. 

In our case, the relatively narrow linewidth of our 
556 nm probe laser (ss 2ir x 10 kHz) allows us to resolve the 
line shape of the excitation spectrum. From the 1 S'o — > 3 Pi 
spectra we extract the Gaussian contribution to the line 
profile, thus measuring the sample temperature. Rabi 
broadening and light shift are also measured and their 
values are compared with the results of a simple model 
describing the 1 Sq — > : Pi trapping transition. We then 
show how time-of-flight temperature measurement can be 
carried out with the same technique. 

We find that the measured temperature is systemati- 
cally above the standard Doppler limit, as observed in 
[221 [23]. We compare our results with the model devel- 
oped in [2T] . We also discuss how this technique could 
give a better insight into the dynamics of trapping and 
cooling of two-electrons atoms. 

In the last section, we present a measurement of the 
cloud temperature and size taken in only 1 ms without 
disturbing the MOT dynamics. This result puts in evi- 
dence that the method discussed in the present paper can 
be regarded as a tool to probe the "instantaneous" cloud 
properties in a non-distructive way. 




Figure 2: Schematic view of the vacuum chamber, (a) 
atomic beam oven; (b) atomic beam probe; (c) differen- 
tial pumping stage; (d) shutter; (e) transverse cooling area 
(not used in the present work); (f) gate valve; (g) Zeeman 
slower; (h) main chamber; (i) ion and non-evaporative 
getter (NEG) pumps; (MOT) MOT beams; (ZS) Zeeman 
slower beam. 
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2 Interpretation and analysis of 
excitation spectra 

Our apparatus is schematically represented in figure [2] In 
the main chamber a 174 Yb magneto-optical trap (MOT) is 
operated on the (blue) 1 Sq — > 1 P\ transition. The 399 nm 
source is a frequency-doubled diode laser stabilized to a 
hollow cathode lamp HI]. The MOT is constituted 
by three retroreflected beams (10 mm diameter, 2.7 mW 
maximum power). About 10 mW are used for the Zee- 
man slower which decelerates thermal atoms from the oven 
(500°C) into the capture range of the MOT. 




Figure 3: Schematic view of the green laser setup. The 
1112 nm IR light is frequency doubled by single passage 
through a ppKTP crystal. Excitation spectroscopy on the 
atomic beam — top left, see also figure part b — pro- 
vides the error signal for locking the 556 nm frequency to 
the 1 S'o — > 3 Pi transition. A double passage Acousto- 
Optic Modulator (AOM2) is used to vary the probe beam 
frequency. 

While the atoms are confined in the blue MOT we excite 
the 1 S'o — » 3 Pi transition using probe light overlapped 
with the MOT beams. Each probe beam has 40 /LtW power 
and 7 mm diameter. The 556 nm light is generated by 
frequency-doubling of a 1112 nm fiber laser in a ppKTP 
crystal. The probe linewidth of 2tt x 10 kHz ~ r p /10 
has been measured by analyzing the transmission of an 
High-Finesse (J- — 5 x 10 4 ) Fabry-Perot resonator (length 
4.74cm and linewidth 2tt x 70 kHz). The position of the 
resonance transition used as reference for the laser lock 
circuit, is determined within 2ir x 200 kHz. This sets the 
limit of our absolut frequency calibration. The green laser 
setup is detailed in figure [31 

As the green probe frequency is scanned we detect 
the 556 nm fluorescence light with a photomultiplier tube 
(PMT). A typical excitation spectrum is shown in figure [4] 



The relevant levels and transitions are sketched in figure 
[5] The three peaks visible in the signal correspond to the 
three \g) — > |po,±i) excitation channels [33]. The width of 
the \g) — ¥ \pq) excitation process (central peak) is deter- 
mined by three independent effects: (i) the bare linewidth 
T p of the 1 S'o —> 3 Pi transition, (ii) the Rabi broadening 
of the 1 S'o level induced by the presence of the 399 nm 
trapping light [27J and (iii) the Doppler broadening due 
to the finite sample temperature. The cloud temperature 
T is obtained from the Gaussian width ad of the Doppler 
contribution using the following relation: 

k B T = -m\la% (1) 

where ks is the Boltzmann constant, m is the atomic mass 
and X p — 556 nm is the probe wavelength. 
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Figure 4: Typical excitation spectrum obtained by prob- 
ing the So — > 3 Pi transition on Yb atoms trapped in 
a 1 < S'o — > 1 Pi MOT and measuring the fluorescence rate. 
The central peak is originated by the \g) — > \po) excitation 
process. The \g) — > \p±i) channels are responsible for the 
two side peaks. See figure [5] for a detailed description of 
the notation used in labeling the relevant energy levels. 

In order to extract the relevant information contained 
in the excitation spectrum we fit the central peak with a 
Voigt function, using the peak position and the Lorentzian 
and Gaussian widths as fitting parameters. The depen- 
dence of the line center on the MOT laser power, gives the 
light shift of the \g) level. The Lorentzian contribution is 
the sum of the natural linewidth and the Rabi broadening 
of the ground state, while the Gaussian width is given by 
the sample temperature. The temperature measurement 
is reliable only if the fitting procedure separates correctly 
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the Gaussian and lorentzian contributions. For this reason 
it is important to properly describe the effect caused by 
the 1 S ^ 1 Pi MOT light onto the 1 S -)• 3 Pi probe and 
compare them with the experimental results. 




Figure 5: (Colors online). Schematic representation of the 
energy levels involved in the 1 Sq — > 3 Pi fluorescence mea- 
surement. \g) stands for the 1 Sq ground state, eo,±i) for 
the 1 P\ multiplet and |po.±i) stand for the m = 0, ±1 
Zeeman sublevels of the 3 Pi state. The \g) — >■ |e) trap- 
ping light (blue continuous arrow) induces a light shift on 
the \g) ground state. The probe beam (green dashed ar- 
row) excites the \g) — > |po,±i) transitions producing the 
characteristic three peak structure visible in the excita- 
tion spectrum. The \p±i) states are Zeeman shifted by 
the quadrupole field used for the MOT. 

We model the \g) — > |e) transition as a four level system 
(taking into account the three Zeeman sublevels |eo,±i)) 
neglecting the presence of the probe beam. The effective 
Hamiltonian is therefore: 
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where At and fit are the trap detuning and Rabi fre- 
quency, respectively. To simplify the notation we intro- 
duce the quantity A f = At + iT t . The Zeeman splitting 
is given by Auj' z . The Rabi frequencies on all three transi- 
tion are assumed to be equal since on average an atom in 
the MOT sees unpolarized light. We find the ground state 
energy E g by diagonalizing this Hamiltonian. The imagi- 
nary part Im(Eg) describes the broadening of the ground 
state induced by laser excitation (Rabi broadening). On 
the other hand Re(E g ) accounts for the energy shift of \g) 



induced by the MOT beams (light shift). The agreement 
between the measured light shift and the result of this 
model is a prerequisite for the validity of the temperature 
measurement. 

For the \g) —t \p±i) transitions (side peaks) similar con- 
siderations are applied. The three effects described above 
contribute to the spectral line broadening in the same 
way. However, for these transitions the presence of the 
MOT magnetic field plays a crucial role. In fact the \p±i) 
states are Zeeman shifted by the trapping quadrupole 
field, which combined with the finite sample size gives rise 
to an additional line broadening. This contribution to the 
transition linewidth reads: 



Li B g P m p t 



(3) 



where [Ib is the Bohr magneton, g p — 3/2 is the Lande fac- 
tor for the — > 3 Pi transition, m p = ±1 is the magnetic 
quantum number, b' is the MOT field gradient and w is the 
cloud size (such that the spatial part of the atomic prob- 

r 2 i 

ability distribution is proportional to exp — ). The 
Zeeman effect, together with gravity, is also responsible 
for the frequency shift of the \g) — > \p±i) lines: due to 
imperfect MOT beam balance and alignment, the atomic 
cloud is formed off center with respect to the quadrupole 
field minimum. The resulting non-zero average magnetic 
field introduces an offset in the resonance frequency equal 
to: 

fiBg P m p b'ro 



Auj 2 = 



(4) 



ro being the center position of the atomic cloud with re- 
spect to the quadrupole field zero. 

3 Experimental measurement of 
line shift and broadenings 

As described in the previous section, the excitation spec- 
trum obtained by probing the So — > 3 Pi transition car- 
ries information about the properties of the atomic cloud 
as well as the characteristics of the trapping potential. In 
the present section we apply this technique to the mea- 
surement of these quantities as a function of the MOT 
saturation parameter St — It/I s ,t- 

During the experiment we keep the trap detuning A t = 
0.86 r t and the magnetic field gradient b' — 36 G/cm fixed. 
While the atoms are trapped in the MOT the probe detun- 
ing A p periodically varies from — 2tt x 40 to +2n x 40 MHz 
in T = 10 ms. For each choice of St we record the 556 nm 
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fluorescence with a PMT for 1 s (i.e. one trace consist of 
100 spectra). Since the AOM efficiency is not uniform over 
the whole scanning range we divide the PMT signal by the 
probe power (monitored by a photodiode during the ac- 
quisition). This operation does not affect the line-shape of 
the spectral feature we are interested in, since the AOM 
efficiency loss is only significant at the borders of the scan 
region. However, it is crucial to have a flat background. 

Once the fluorescence has been acquired and processed 
we extract from the full waveform the individual spectra 
and fit the central peak with a Voigt profile. We average 
the fit results and assign to each parameter an error bar 
equal to the calculated standard deviation. 

3.1 Ground state light shift 

We first focus our attention on the ground state light shift. 
The experimental results are shown in figure [5] where the 
\g) — >• |po) peak shift Au> g is plotted as a function of the 
MOT saturation parameter. A direct comparison of the 
experimental data with the model presented in section [5] 
is valid only assuming that the horizontal axis calibration 
is not affected by systematic errors and the resonance con- 
dition is properly determined. Since we can not know a 
priori if these conditions are fulfilled we proceed in a dif- 
ferent way. 

We diagonalize the Hamiltonian of eq. [2] using At = 
0.86 r t and Alo' z = 2tt x 12.4 ± 0.4 MHz. This last 
value is calculated from the side peak splitting Au z = 
2ir x 18.6 ± 0.6 MHz observed in the experiment by using 
Auj' z = Au z gt/g p , being g t = 1 and g p = 3/2 the Lande 
factors for the trap and probe transition respectively. We 
fit the experimental data for Aw 9 with Ke(E g )/h derived 
from the diagonalization. As fitting parameters we use the 
calibration factor a (i.e. the saturation parameter used in 
H is ast) and the offset Aw g o- We obtain a — 0.67 ±0.07 
and Awgo = 2-7T x 410 ± 80 kHz. The calculated behaviour 
is represented in figure [6] as a continuous line. 

The fact that Aw g o is different from zero indicates that 
the 556 nm laser source may be locked with non-zero de- 
tuning with respect to the atomic resonance. However 
it also has to be considered that during the measurement 
process the Zeeman slower beam is constantly present thus 
contributing to the ground state light shift. With only 
this measurement we cannot distinguish between the two 
effects. 

Concerning the value of a it is important to notice 
that the saturation parameter calibration assumes that 
the atomic cloud is formed at the intersection of the six 
trapping beams, considered to be perfectly Gaussian. It is 
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Figure 6: Position of the \g) — > \po) peak relative to the 
resonance condition as a function of blue MOT saturation 
parameter s t for A t = 0.86 T t . The continuous line is 
obtained from the four-level model presented in section 
[21 The fitted offset and calibration factor are Au g0 = 
2tt x 410 ± 80 kHz and a = 0.67 ± 0.07 respectively (see 
text for details). 

realistic to assume that the actual intensity at the MOT 
position is lower than the estimated value. Therefore the 
measured mismatch of about 30% can be considered to be 
in reasonable agreement with our expectation. 

3.2 Mapping the trap laser intensity: an 
example 

The observed value of a gives us information about the 
actual light intensity at the MOT position. This suggests 
that the trapped atoms can be used as a probe to explore 
the spatial distribution of the light field for a given MOT 
alignment. 

The measurement is performed as follows. First we load 
the 1 S -> 1 P 1 MOT and we probe the 1 S -> 3 Pi tran- 
sition recording the excitation spectrum. We then shift 
the quadrupole field zero by applying a uniform mag- 
netic field and we perform a new fluorescence measure- 
ment. The sequence is repeated for different trap positions 
while the laser beam alignment remains unchanged. Two 
pairs of Helmholtz coils generate a uniform magnetic field 
along the vertical (x) and Zeeman slower (y) directions of 
8.6 G/A and 4.1 G/A respectively. For the typical gradi- 
ent of 36 G/cm (18 G/cm along y) used in the experiment 
the corresponding position shifts are 0.28 mm/G along x 
and 0.56 mm/G in the y direction. 

During the mapping sequence we vary both coil cur- 
rents from —0.5 to 0.5 A in steps of 0.1 A. For each cur- 
rent configuration we acquire 10 excitation spectra and we 
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Figure 7: Position-dependent ground state light shift 
ameasured through variation of compensation coil cur- 
rents. The frequency shift is referred to its value in ab- 
sence of compensation field. Laser intensity and detuning 
are maintained fixed at s t — 0.5 and A t = 0.86 T t . 

analyze them as described above. In figure [7] we plot the 
frequency offset of the \g) — > \po) peak as a function of 
the offset fields. Here Aw 9 is referenced to its value in 
absence of external fields (i.e. Aui g = for B x ,B y = 0). 
During the experimental sequence the MOT light inten- 
sity was kept constant at St — 0.5. Therefore the observed 
light shift variation is only due to the displacement of the 
MOT across the intensity distribution of the trapping light 
field. 

We observe that the light shift increases with B x (i.e. as 
the atoms are moved downwards). Therefore the light field 
maximum is located below the cloud position for B x = 0. 
Since we aligned the MOT for B x = B y = maximizing 
the number of atoms, we conclude that for the optimum 
alignment the cloud is formed above the intersection of the 
trapping beam. We also observed that the optimum MOT 
loading is achieved when the Zeeman slower beam passes 
below the atom position. An interpretation of this result 
would require a theoretical treatment which is beyond the 
purpose of this paper. 

The experiment presented here is an example of how 
the technique discussed in this paper can be used to mea- 
sure the properties of the trapping potential. For this 
particular measurement we observe the \g) — > |po) peak 
shift and we relate it to the local trap field intensity. In 
the following we will show how a proper determination 
of the Lorentzian and Gaussian contributions to the peak 
linewidth, can provide a reliable measurement of the sam- 
ple temperature. We think that this technique could be 
helpful for unveiling the actual cloud dynamics. 



3.3 Rabi broadening 

The presence of 399 nm MOT light also modifies the 
linewidth of the ground state \g) (Rabi broadening). As 
a result the Lorentzian contribution to the central peak 
width (T gp ) increases as a function of the trap saturation 
parameter St- In figure [8] we plot the experimental values 
obtained from the data set of figure [6] The comparison 
with the four-level model is carried out following the same 
procedure as for the ground state light shift analysis. In 
this case the Hamiltonian of eq. [H is diagonalized assum- 
ing that a — 0.67. The only free parameter is therefore 
r gp oj i-G- the value of T gp in the limit of zero intensity. 
A minimization process similar to the one applied before 
leads to T gp0 = 2vr x 800 ± 60 kHz. The fit is plotted in 
figure [5] as a continuous line. 
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Figure 8: Central peak Lorentzian width T gp as a function 
of MOT saturation parameter St- The continuous line 
is obtained diagonalizing the Hamiltonian of eq. [2] with 
a = 0.67. The minimization process gives an offset T gp o = 
2tt x 800 ± 60 kHz (see text for details). 

As discussed in section [2] the Lorentzian width of the 
\g) — > \po) peak is also determined by the bare atomic 
transition linewidth. Therefore for s t = we expect to 
have T gp Q — T p , where T p — T p ^l + s p . From the probe 
beams size and power we have s p ~ 9, therefore r p m 2ir x 
580 kHz. The difference between T gp o and T p is attributed 
to the presence of the Zeeman slower beam during the 
measurement and to the linewidth of our probe laser (2ir x 
10 kHz). For st > 0.25 the experimental data show good 
agreement with the fit. The deviation observed for st < 
0.25 is due to the fact that the Voigt fit is not accurate 
when the Lorentzian contribution to the line profile is too 
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small. 

3.4 Cloud temperature 

The analysis described above shows that both light shift 
and Rabi broadening are consistently taken into account 
and extracted from the fluorescence profiles. The residual 
Gaussian width can now be related to the cloud tempera- 
ture via eq. [T] Experimental data are plotted in figure [9j 
The continous line represent the standard Doppler limit 
Td = 695 /iK which is independent of the trap light in- 
tensity. The dot-dashed line is derived from the model 
presented in the paper by Choi et al. (eq. 13 of [2"T]). 
The dashed line is obtained in the same way but replacing 
the experimental value of St with ast being a = 0.67 the 
calibration factor derived from our previous analysis. 
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Figure 9: Cloud temperature as a function of St. The 
continuous curve is the standard Doppler limit, the dot- 
dashed and dashed line are obtained from the theoretical 
model discussed in [21] calculated with a = 1 and a = 0.67 
respectively (see text for details). 

We observe that the measured temperature is higher 
than the standard Doppler limit as predicted by the model 
of reference [21]. However experimental data are system- 
atically above the theoretical curve and when the calibra- 
tion factor is applied the deviation between experiment 
and theory increases. The lack of quantitative agreement 
with the experiment could be related to the presence of 
other effects not included in the model. In particular, as 
shown in [3D], spatial intensity fluctuations of the cool- 
ing beams are responsible for extra- heating. From this 
perspective the knowledge about the actual trap light in- 



tensity at the cloud position observed in section I3T21 could 
be a relevant ingredient for a proper theoretical treatment 
of the cooling dynamics of two-electrons atoms. 

4 Time resolved measurement 

As discussed in section[3]the presence of the Zeeman slower 
beam during the measurement process and the detuning 
of the probe laser with respect to the 1 5'o — > 3 Pi transition 
influence the determination of the ground state light shift. 
In order to avoid these effects we proceed as follows. 

As a first step we perform a fluorescence measurement 
on trapped atoms as described in the previous section. We 
then switch off the blue light by means of an AOM. The 
atoms decay in a time r t = l/r t = 5.5 ns to the I 1 So) 
state. Since the ground level has angular momentum, 
the MOT gradient does not influence the atomic motion. 
Finally we perform a second fluorescence measurement on 
free falling atoms. 

The absence of confinement during the last stage forces 
us to realize the measurement faster than the time needed 
for the atoms to leave the probe region. For this reason 
we keep the frequency span at 2ir x 80 MHz as in section 
[3] and reduce the scan period to 1 ms. Note that during 
this time the cloud's center-of-mass moves by about 5 /Ltm. 
With the resulting scan rate of R — 2ir x 80MHz/ms the 
probe frequency stays on resonance for a time equal to 
t s = T p /R = 2.3 /is. This value has to be compared to the 
atomic decay time r p = 870 ns. The fact that t s « 2.6t p 
ensures that the fluorescence signal is not affected by non- 
adiabaticity. 

The difference between the \g) — > \po) peak positions 
with and without trapping light is only related to the 
ground state light shift caused by the MOT beams. The 
experimental results are shown in figure [10] (bottom) to- 
gether with the model comparison (see section [3]) . The 
minimization process gives us a — 0.55 ±0.04 and Aw 5 o = 
27r x 30 ± 60 kHz. The calibration factor a is in agreement 
with the previously measured value. The offset Aw g o is 
compatible with zero as expected. 

From the same data set we obtain the ground state Rabi 
broadening as described in section [3] For this purpose we 
use the spectra acquired in the presence of 399 nm light. 
Experiment and model are shown in figure [TD] (top). As in 
section [3] a is kept fixed at the previously obtained value 
of 0.55 while for T gp0 we find 2tt x 800 ± 30 kHz. 

In figure [11] we plot the cloud temperature calculated 
from eq. [T] using the fluorescence signal measured before 
(empty circles) and after (full squares) switching off the 
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Figure 10: Rabi broadening (top) and light shift (bottom) 
as a function of st (vertical axis unit is MHz). The ex- 
perimental data are obtained comparing spectroscopy in 
the presence and absence of 399 nm trapping light. The 
continuous lines are obtained from the four-level model of 
section [2] The trap detuning is A t = 0.86 T t . The mag- 
netic field gradient is 36 G /cm. 



Figure 11: Cloud temperature measured with (empty cir- 
cles) and without (full squares) MOT light. For the empty 
circles the experimental errors are smaller that the symbol 
size. The data of figure [TT] are also shown (empty trian- 
gles). The continous line is the standard Doppler limit. 
The prediction given by the model of [21] is plotted as a 
dash-dotted line (dashed with calibration correction). 



blue MOT. In the latter case the spectra are fitted with a 
fixed Lorentzian width equal to the value of T gp o obtained 
from our previous analysis. The two sets are compatible 
showing that the Voigt fit correctly identifies the Gaussian 
component independently of the presence of the trapping 
light. As for figure [9] we plot the standard Doppler limit 
(continuous) and the theoretical curve after reference [21] 
calculated for a = 1 (dot-dashed) and a = 0.55 (dashed 
line). 

As detailed above, each excitation spectrum is acquired 
in T = 1 ms. This value is 10 times smaller than the 
one used for the experiments described in section [3] How- 
ever, if the scan is adiabatic, the temperature measure- 
ment does not depend on the particular choice of T. In 
order to verify this assumption, we reproduce in figure [TT] 
the experimental points previously presented in figure [9] 
(empty triangles). We note that, whitin the experimental 
error, the three data sets. 

Due to the relatively short scan period T used in the 
experiment, we regard the information contained in each 
excitation spectrum as "instantaneous" . In other words 
we are performing a time resolved measurement of the 
atomic cloud properties. This suggest a different way of 
comparing the temperature obtained using the technique 



of [T] [2] with the result of a more standard time-of-flight 
(TOF) approach. The procedure we use is described as 
follows. 

We first load the 399 nm MOT and then we switch off 
the trapping beams. After a variable time t we perform a 
single probe frequency scan in 1 ms. We fit the three peaks 
present in the fluorescence signal with Voigt profiles. We 
then calculate the difference between the Gaussian widths 
of the central and side peaks. In this way we obtain the 
contribution a z to the side peak broadening due to the 
presence of the MOT magnetic field gradient. Finally a z 
can be related to the cloud size via equation [3] This pro- 
cedure resembles the standard TOF technique. The only 
difference between the two approaches resides in the way 
used to determine the sample size: in one case by ab- 
sorption/fluorescence imaging (direct), in the other one 
by excitation spectroscopy (indirect). 

In figure [T2l we plot the cloud size w as a function of the 
free fall time t. As it is clearly visible in the graph w is con- 
stant for t > 2.5 ms. This effect can be easily understood 
considering that as t increases the cloud becames bigger 
and eventually its image on the PMT exceedes the detec- 
tor aperture. When this occurs only part of the atomic 
distribution contributes to the observed fluorescence sig- 
nal. In figure [TT] we place a dot-dashed horizontal line at 
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a position equal to one half of the PMT aperture radius. 
When the actual cloud size is equal to this value, 5% of 
the atomic density ditribution (defined as in section [2]) is 
not imaged onto the detector. 
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Figure 12: Cloud size w as a function of TOF t. The hor- 
izontal dot-dashed line represent the detector aperture. 
The continuous line is the obtained fitting the experimen- 
tal data with the function of eq. [5j 

We fit the data points for t < 2.5 ms with the expected 
expression for free ballistic expansion: 

w(t) = W (0W1 + — i 2 , (5) 

V 771 

where m is the atomic mass and w(0) is the initial cloud 
size. From the fit we get T = 5.5±1.2mK. This result has 
to be compared with the correspondent value determined 
by a fully spectroscopic method. For this purposes we 
analyze the spectra acquired before releasing the trap fol- 
lowing the same procedure as in section [3] This procedure 
leads to T = 4.91 ± 0.06 mK which is in good agreement 
with the TOF value. 

5 Conclusions 

We extended the method proposed in [TJ [2] for online 
monitoring of the temperature of cold atoms. Probing 
the 1 5o — ) 3 Pi transition of 174 Yb atoms confined in a 
magneto-optical trap we measured the ground state light 
shift and broadening induced by the MOT light field. We 
then fit the experimental data with a simple model de- 
scribing the 1 Sq — > 1 Pi transition. As a results we cal- 



ibrated the MOT local intensity. This allowed us to ex- 
tract the Gaussian component of the excitation spectra 
and measure the sample temperature avoiding the limita- 
tions encountered by the authors of [2J. 

We applied a similar procedure to free-falling atoms. We 
compared the measurements performed with and without 
confining potential thus obtaining the ground state light 
shift in absence of spurious effects (like additional light 
shifts and laser locking offset). We showed that the tem- 
perature measurement is not affected by 399 nm light indi- 
cating that the Gaussian contribution to the line is prop- 
erly determined. Finally we measure the sample temper- 
ature from free cloud expansion. In contrast with tradi- 
tional time-of- flight measurement we used the \g) — > \p±i) 
linewidth for determining the cloud size. We showed that 
the results obtained with the two method are in agree- 
ment. 

The measured temperature of our Yb MOT is consis- 
tently above the Doppler limit. This result agrees with the 
observations performed by other groups [22j [23] . We con- 
firmed that the theoretical model presented in [2T] qualita- 
tively describes the dependence of the cloud temperature 
on the trap light intensity. However our data do not show 
a quantitative agreement with the theory. We measured 
the actual MOT light intensity distribution at the cloud 
position. We pointed out that this could account for the 
observed disagreement. 

Providing that the adiabaticity condition is satified, the 
acquisition time T can be reduced without affecting the 
information carried by the excitation spectrum. The value 
of T = 1 ms used in the experiment, does not have to be 
regarded as a fundamental limit. Keeping the scan rate R 
fixed, T can be lowered by choosing a different frequency 
span. Thanks to this fact the cloud temperature and size 
could be monitor on a fast time scale. 

In conclusion the technique presented in this paper 
opens the possibility to monitor online and in "real time" 
the dynamics of a cold atomic cloud. This could be a 
key tool to understand the peculiar thermodynamics of 
magneto-optical trapping of two-eletrons atoms. 
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